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! .0  SUMM  ARY 


I Ins  document  summarizes  the  symmetric  flutter  analyses  conducted  on  the 

;I"  l1,U)V:lblc  hor,/ont:i1  t;"'  <)r  thL>  B‘>cing-designed  SSI  The  analytical  work  was  completed 
undci  contract  to  the  Department  ol  I ransportation  as  a follow-on  technology  program 
alter  SSI  program  cancellation.  The  horizontal  stabilizer  was  analyzed  lor  flutter 
instabilities  including  the  interaction  effects  of  the  fuselage  and  wing  structure,  the  airplane 
contiol  systems,  and  the  elevator  gearing  ratio. 

Results  of  the  analyses  are  given  in  sections  4.4  and  5.3.  Rased  on  these  results  it  is 
concluded  that  the  SSI  empennage  provided  adequate  margins  of  safety  for  flutter  lor  all 
flight  conditions  (including  failure  conditions)  except  when  the  stabilizer  was  blown  hack 
due  to  high  aerodynamic  loading.  It  is  also  concluded  that  the  flutter  speed  ol  an 

all-movable  horizontal  stabilizer  with  geared  elevator  is  relatively  insensitive  to  elevator 
gearing  ratio. 

I he  analyses  reported  in  this  document  incorporated  linear  representation  of  some 
nonlinear  and  transient  control  system  conditions.  Although  this  is  considered  adequate  to 
account  tor  lust-order  ellects  on  flutter  speeds,  an  extension  of  these  analyses  to  define  the 
nonhnear  ellects  is  recommended.  Such  studies  should  employ  analog-digital  (hybrid) 
computer  techniques  lor  the  evaluation  of  transient  actuator  conditions. 


2.0  INTRODUCTION 


I ollowing  cancellation  «l  tlie  SSI  Program  in  1 ‘>7 1 , a program  was  initiated  to 
complete  several  Hatter  analyses  of  the  liori/  mtal  tail  to  preserve  technical  information  lor 
design  support  and  development  of  future  airplanes  of  this  general  type.  I he  program 

consisted  ol  completion  ol  the  symmetric  flutter  analyses  of  the  empennage  structure  ol 
the  SST. 

the  horizontal-tail  design  incorporated  .some  unique  features  of  interest  and 
importance  to  future  designs  of  this  type.  These  included  I I > an  all-movable  slab  stabilizer 
with  geared  elevators.  (2)  a quadruple-redundant  hydraulic  actuation  system  for  the 
horizontal  stabilizer,  (3)  reliance  on  actuator  stiffness  rather  than  surface  mass  balance  for 
the  suppression  ol  Mutter,  and  (4)  the  possibility  that  the  ..erodynamic  hinge  moments 
could  become  large  enough  to  stall,  and  even  back-drive,  the  stabilizer  actuators  with  a 
resultant  substantial  reduction  in  actuator  stiffnesses.  I urthermore,  it  was  necessary  to 

account  lor  the  response  of  the  active  longitudinal  control  system  in  the  stabilizer  Mutter 
analysis 

Initial  empennage  flutter  studies  on  the  SST  represented  the  system  as  empennage 
surl aces  attached  to  the  alt  body,  which  was  cantilevered  from  the  wing  rear  spar  location. 
Although  the  stabilizer-elevator  structural  model  was  quite  detailed,  being  derived  through 
enrol  I in  ite-elenient  representation  of  the  structure,  the  fuselage  itself  was  represented  as  a 
x-ant  and  the  attachment  bulkhead  structure  was  in  a preliminary  stage  of  design. 

While  this  analy  tieal  model  was  considered  adequate  for  preliminary  design  evaluation 
H did  not  lully  represent  the  effects  on  flutter  of  the  wing/body  modes  in  the  frequency 
region  ol  stabilizer  Mutter,  f lutter  speed  boundaries  associated  with  initial  strength  design 
showed  significant  deficiencies,  as  illustrated  in  figure  I.  Addition  of  stiffness  to  the 
hydraulic  actuators,  supporting  structure,  and  stabilizer  carry-through  structure  raised  the 
speeds  to  acceptable  levels  when  the  actuators  were  in  a holding  condition,  as  shown  in 
ligure  I,  but  actuator  failure  conditions  or  the  transient  conditions  of  stabilizer  under 
command,  stalled  actuators,  or  (down-back  actuators  still  produced  Mutter  speed  defi- 
ciencies. Analog  Mutter  studies  were  made  that  carefully  simulated  the  non-linear  effects  of 
external  loading  and  control  inputs  on  dynamic  stiffness  of  the  actuator  system,  using 
quasi-steady  aerodynamic  forces  that  produced  Mutter  speeds  closely  approximating  the 
results  ol  the  linear  analyses.  The  nonlinear  behavior  of  the  actuator  system  did  not 
sigml icantly  alter  the  Mutter  speed  boundaries  associated  with  zero  structural  damping. 

Incorporation  of  a more  detailed  representation  of  the  body  and  inclusion  ol  the  wing 
structural  and  aerodynamic  effects  in  the  linear  analysis  were  underway  at  the  time  of  the 
SSI  program  cancellation  in  the  spring  of  1971,  Results  of  the  completion  of  these  studies 
arc  reported  m this  document.  An  additional  study  related  to  elevator  gearing  effects  on  the 

Mutter  characteristics  ol  the  all-movable  stabilizer  was  conducted  under  the  contract  and  the 
results  are  presented 
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3.0  ANALYSIS  FLOW 


Horizontal-stabilizer  flutter  analyses  using  free-free  airplane  modes  were  conducted 
with  a stillness  matrix  reduced  to  1 94  degrees  ol'  freedom  representing  I tic  symmetric 
half-airplane.  This  stiffness  matrix  was  combined  witli  mass  data  representing  three 
gross-weight  conditions  to  determine  normal  modes  and  frequencies  for  those  cases.  All 
horizontal-stabilizer  actuators  were  initially  represented  as  rigid  links. 

Previous  studies  of  the  analytical  model  employing  a cantilevered  aft  body/horizontal 
tail  indicated  that  a minimum  of  six  stabilizer  modes  were  required  for  an  accurate  flutter 
solution.  Approximately  40  airplane  modes  (depending  on  airplane  gross  weight)  were 
required  in  order  to  cover  the  frequency  range  of  the  first  six  stabilizer  modes. 

Inspection  of  the  calculated  mode  shapes  indicated  that  a few  modes,  particularly  at 
the  higher  frequencies  and  heavier  weight  conditions,  were  not  accurately  representing  the 
structure;  i.e„  individual  mass  resonances  (“tcntpoles”)  appeared  in  local  areas  ol  the  wing, 
l liese  “tcntpoles”  resulted  from  the  necessity  of  distributing  the  mass  data  over  a limited 
number  of  nodes. 

Alternate  mass  redistributions  over  limited  areas  (so  that  the  airplane  eg  was  not 
changed)  proved  to  be  beneficial  in  eliminating  the  “tcntpoles”  in  a number  of  cases.  Those 
modes  associated  with  the  remaining  “tcntpoles”  were  dropped  from  the  analyses  as  being 
invalid  and  not  representative  of  the  actual  structure. 

After  a representative  set  of  mode  shapes  and  frequencies  was  obtained  for  each  weight 
condition,  generalized  aerodynamic  forces  were  determined  at  one  supersonic  and  one 
subsonic  Mach  number  for  each  weight  condition.  Three-dimensional,  compressible 
aerodynamic  theories  (ref.  1 and  2)  were  used  at  M = 0.9  and  1.65.  Flutter  solutions 
including  V-g  and  V-co  plots  were  then  obtained  for  these  conditions.  The  number  of  modal 
degrees  of  freedom  in  these  solutions  ranged  from  34  for  the  lightweight  airplane  to  40  for 
the  heavyweight  airplane. 

These  analyses  served  to  determine  which  of  the  three  weight  conditions  analyzed  was 
critical  from  a Rutter  speed  standpoint.  From  these  results  it  was  decided  to  continue  the 
Rutter  studies  for  only  one  weight  condition,  namely  the  41 7,400-lb  gross-weight  airplane. 

Studies  were  conducted  to  determine  the  necessary  modes  to  be  included  in  the 
remaining  analyses  (deleting  modes  that  had  no  influence  on  the  critical  Rutter  condition). 
The  primary  objective  was  to  reduce  computer  time  requirements  when  generalized 
aerodynamic  forces  were  calculated  for  other  Mach  numbers.  The  major  requirement  of  the 
retention  studies  was  to  minimize  any  change  in  the  stabilizer  flutter  speeds.  The  Rutter 
speed  for  the  basic  wing  was  more  sensitive  to  modal  deletions  than  the  stabilizer  mode, 
with  the  result  that  the  final  set  of  retained  modes  produced  a moderate  change  in  wing 
flutter  speeds  in  some  cases,  although  the  basic  wing  Rutter  mode  and  its  coupling 
characteristics  were  essentially  unchanged  in  all  analyses. 
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Results  of  t lie  modal  retention  studies  showed  that  l<S  modes  ot  the  airplane  were 
required  supersonically  and  26  modes  subsonically  to  accurately  assess  the  horizontal- 
stabilizer  flutter  characteristics.  Additional  flutter  analyses  were  then  conducted  at  M - 0.6, 
I 2,  and  2 1. 


To  assess  the  control-system  effects  on  horizontal-stabilizer  flutter,  the  following  steps 
were  taken  The  aerodynamic  stiffness  and  damping  coefficients  were  determined  lor  several 
reduced  frequencies  in  the  usual  manner.  The  flutter  equations  were  rewritten  in  terms  ol 
the  Laplace  transform.  S,  and  merged  with  the  additional  actuator-mode  and  control 
equations.  The  merged  equations  were  then  solved.  The  only  valid  roots  ol  the  solution 
occur  when  S = ito.  and  the  locus  of  these  roots  (as  iterations  were  made  on  the  reduced 
frequencies)  determined  the  neutral  stability  (flutter  point).  Hie  details  of  this  formulation 
are  given  in  appendix  A.  By  changing  the  transfer  function  (stiffness)  of  the  horizontal- 
stabilizer  actuators,  various  actuator  operating  conditions  such  as  holding,  stalled, 
blown-back,  etc.,  were  readily  assessed. 


I valuation  of  the  failed-aetuator  conditions  was  accomplished  in  the  same  manner  as 
for  the  non  failed  conditions.  I lie  Iree-tree  airplane  stillness  matrix  was  modified  by 
removing  the  link  representing  the  outboard  actuators.  New  modes  were  calculated  and  the 
analyses  proceeded  in  a manner  similar  to  that  for  the  nonfailure  configuration. 


4.0  FREE-FREE  AIRPLANE  ANALYSIS 


4.1  DESCRIPTION  OF  STRUCTURE 

The  general  configuration  of  the  USA  supersonic  transport  is  shown  in  figure  2.  The 
stiffness  matrix  used  for  the  empennage  flutter  studies  was  determined  from  a ^ 

representation  of  the  symmetric  half-airplane,  using  beam  and  plate  elements.  urc  3 
shows  the  horizontal-stabilizer  representation.  Initially,  over  2000  nodes  were  used 
represent  the  airplane  structure.  These  nodes  were  chosen  to  provide  the  greatest  deta 
possible  in  representing  the  empennage  structure;  to  avoid  excessive  computational 
requirements  some  simplifications  were  required  in  representing  the  forward  fuselage  and 
wing  structure.  The  basic  dynamic  characteristics  of  the  wing  and  lorward  fuselage  were  not 
significantly  altered  by  these  simplifications.  The  representation  employed  194  nodes 
distributed  as  follows; body,  19;  wing,  46;  engine  nacelles.  12;  horizontal  stabilizer,  68.  fin 

and  ventral  surface,  49. 

The  horizontal-stabilizer  structure  was  of  skin-spar-rib  construction  in  the  mam  box 
area  with  full  depth  honeycomb  at  the  leading  edge,  tip,  and  elevator.  Figure  4 shows  the 
stabilizer  structural  configuration.  Also  shown  in  figure  4 is  a schematic  representation  ol 
the  elevator  gearing  design. 

4.2  DESCRIPTION  OF  CONTROL  SYSTEM 

The  complete  longitudinal  flight  control  system  is  shown  schematically  in  figure  5 
(reproduced  from  reference  3,  where  it  is  explained  in  more  detail).  Commands  were 
transmitted  from  the  pilot’s  controls  by  mechanical  linkage  and  cables  to  the  master 
servo  was  connected  mechanically  to  the  control  valves  on  the  four  stabilizer  actuators. 
Pitch  commands  were  also  transmitted  electrically  via  four-channel  USAS  Hardened 
Stability  Augmentation  System)  and  ECSS  (Electric  Command  and  Stability  System)  paths 
to  the  four  electric  command  (EC)  servos.  The  output  of  the  EC  servos  was  mechanically 
summed  with  the  output  of  the  master  servo  at  each  surface  actuator.  Multiple  load  pa n 
were  provided  from  the  master  servo  to  the  four  stabilizer  actuators.  Each  o the  fou 
actuators  was  powered  by  a separate,  independent  hydraulic  system.  The  flight  control 
system  used  in  the  flutter  analyses  is  shown  in  the  simplified  block  diagram  of  lljure  6. 

The  SST  all-movable  horizontal  stabilizer,  chosen  for  its  high  degree  of  control 
effectiveness,  was  constrained  against  rotation  solely  by  the  stiffness  ol _ the  actuation 
system.  Figure  7 shows  a typical  actuator  installation.  The  stiffness  of  the  full-scale 
actuation  system  was  determined  in  an  experimental  program  and  reported  on  in  the  Phase  1 
SST  Technology  Follow-On  program  (ref.  3).  The  results  of  these  tests  were  incorporated  in 
the  flutter  analyses  reported  herein.  Figures  8 and  9 show  the  magnitude  of  stif  fness  and  the 
phase  of  actuator  reaction  force  relative  to  actuator  deflection  when  an  oscillatory  load  is 
applied  to  the  output  end  of  a single  horizontal-stabilizer  actuator  under  normal  holding 

conditions. 


The  SSI  liorizontiil-stabilizer  design  was  swell  that  the  surface  liinge  moments  were 
large  enough  to  stall  the  actuators  under  certain  extreme  boundary  emergency  (light 
conditions,  resulting  in  a significant  decrease  in  actuator  stiffness.  The  softness  decrease 
results  from  the  main  servo  valve  being  wide  open,  thus  allowing  relatively  long  hydraulic 
lines  to  become  the  dominant  elastic  element  in  the  system.  The  actuator  stiffness  under 
stalled  conditions,  as  determined  in  reference  3,  is  shown  in  figure  10.  The  phase  shift  is 
shown  in  figure  II. 

A possibility  existed  that  an  extreme  airplane  maneuver  in  atmospheric  turbulence 
could  produce  forces  that  would  exceed  actuator  capability.  Under  this  condition,  the 
actuators  would  be  driven  back  from  their  commanded  position,  resulting  in  a significant 
decrease  in  actuation  system  stiffness.  Both  the  main  servo  valve  and  the  relief/bypass  valves 
are  open  under  blowback  conditions.  The  actuator  stiffness  and  phase  shift  under  blowback 
conditions,  as  determined  by  the  reference-3  test,  are  shown  by  the  solid  lines  in  figures  12 
and  13.  This  drastically  reduced  stiffness  was  of  major  concern  during  the  design  of  the  SST 
since  it  produced  a statically  divergent  mode  at  extremely  low  speeds.  Immediately  prior  to 
the  SST  program  cancellation  one  alternate  design,  involving  a time-delayed  bypass  valve 
operation,  was  presented  for  flutter  evaluation.  Its  transfer  function  is  shown  by  the  dashed 
lines  in  figures  12  and  13.  Results  of  the  actuator  stiffness  tests  of  reference  3 are  reflected 
in  the  flutter  studies  of  this  document,  with  the  exception  of  the  proposed  alternate  design, 
which  was  not  checked  during  the  full-scale  hardware  test  of  Phase  1. 

For  high  command  rates  the  stabilizer  main  servo  valve  is  nearly  wide  open,  (i.e.,  is  at 
approximately  the  valve  opening  for  the  stalled  actuator  position).  At  low  command  rates 
the  main  servo  valve  is  nearly  closed,  (i.e.,  is  at  approximately  the  valve  opening  for  the 
holding  actuator  position).  Since  the  under-command  stiffness  is  bounded  by  two  stabilizer 
operating  conditions  (stalled  and  holding)  that  were  analyzed  fully,  no  analyses  were 
conducted  for  under-conunand  conditions. 

4.3  DESCRIPTION  OF  WEIGHT  CONDITIONS 

Paneled  weights  data  for  three  configurations  were  employed  initially  representing 
airplane  gross  weights  of  417,400,  529,900  and  601,9001b.  The  weight  conditions  also 
represent  a range  of  center-of-gravity  locations  approaching  fore  and  aft  extremes  of 
expected  flight  conditions  for  the  airplane.  Table  1 presents  weight  condition  detail. 

Of  particular  concern  in  the  empennage  flutter  studies  was  the  fact  that  the  SST  design 
incorporated  an  aft  body  fuel  tank  (designated  8A).  Previous  empennage  flutter  studies  on 
the  SST  had  indicated  that  the  fuel  level  in  this  tank  could  have  a significant  effect  on  the 

empennage  flutter  speeds.  The  location  of  this  tank  in  relation  to  the  empennage  structure  is 
shown  in  figure  4. 

4.4  ANALYSIS  RESULTS 

4.4.1  Stabilizer  Holding  Conditions 

I he  three  airplane  gross  weight  conditions  were  initially  analyzed  for  flutter  at  M = 0 9 
and  1.65  with  actuators  considered  to  be  rigid  links  and  without  control  system  effects 


TABLE  1. -AIRPLANE  WEIGHT  CONDITION  DESCRIPTIONS 


included.  Additional  analyses  were  conducted  at  M = 0.6,  1.20  and  2.1  to  complete  the 
curves  ol  figure  14.  Similar  analyses  were  run  lor  all  three  weight  conditions  with  the 
actuator  dynamic  properties  of  figures  8 and  9 included.  The  results  are  shown  in  figure  1 5. 
I he  llutter  speeds  are  somewhat  lower  than  for  the  rigid-link  cases  but  have  the  same 
general  trend.  The  llutter  speeds  for  the  empennage  as  a part  of  the  lull  airplane  are 
substantially  higher  than  for  the  cantilevered  body  representation  of  figure  1 ; as  a 
consequence,  lor  the  stabilizer  holding  condition  the  flutter  margins  art'  ample.  This  is 
believed  to  be  partly  the  result  of  the  use  of  updated  versions  ol  the  aerodynamics  ol 
references  1 and  2 and  partly  the  icsult  of  the  use  of  a more  detailed  mathematical  model  of 
the  luselagc  structure,  including  the  interface  between  horizontal  hiil  and  body.  Flutter 
studies  ot  elevator  gear-ratio  effects,  section  'TO,  employed  a cantilevered  body  representa- 
tion that  included  the  detailed  mathematical  model  of  the  aft  body  and  stabilizer 
attachment.  These  results  agreed  very  well  with  the  full-span  data,  indicating  that  the 
wing/torebody  did  not  change  the  stabilizer  flutter  speeds  significantly.  The  significant 
change  tn  flutter  speed  from  the  earlier  studies  in  spite  of  relatively  minor  changes  in 
structural  representation  is  believed  to  be  due  to  inaeeuraeies  in  unsteady  aerodynamic 
theory.  This  emphasizes  the  need  for  experimental  verification  of  analytical  results  for 
configurations  of  this  type. 

4.4.2  Selection  of  Critical  Weight  Condition 

At  M = 0.9,  the  lightweight  airplane  condition  (417,400  lb)  produced  the  lowest  flutter 
speeds,  as  shown  in  the  upper  graph  of  figure  16.  However,  when  considering  the  llutter 
margin  available  over  1.15  V^,  shown  by  the  open  symbols  of  the  lower  graph,  the 
intermediate  gross  weight  configuration  (529,900  lb)  was  critical  by  a small  amount.  The 
flutter  characteristics  and  modal  coupling  of  the  417,400-  and  529,900-lb  conditions  were 
very  similar. 

As  shown  in  figure  15  the  lightweight  and  heavyweight  configurations  had  virtually 
identical  flutter  speeds  at  M = 1.65.  Thus,  at  supersonic  speeds,  the  601,900-lb 
configuration  was  somewhat  more  critical  than  the  417,400-lb  weight  condition  because  of 
the  more  stringent  flight  placard.  However,  the  overall  similarity  in  the  characteristics  of 
flutter  of  the  horizontal  tail  for  the  range  of  gross  weights  studied  led  to  the  conclusion  that 
one  weight  condition  was  probably  sufficient  for  further  flutter  studies  of  control  system 
variations. 

Previous  studies  on  an  analytical  model  with  a cantilevered  aft  body  and  horizontal  tail 
indicated  that  the  lowest  stabilizer  flutter  speeds  occurred  when  the  aft  fuel  tank  8A  (fig.  4) 
was  partially  full.  The  data  shown  in  figure  17  for  M = 0.9  for  the  full  span  ease  was  in 
agreement  with  these  previous  results.  This  led  to  the  selection  of  the  417,400-lb 
configuration  for  further  studies. 

4.4.3  Actuator  Failure  and  Extreme  Loading  Conditions 

The  potential  reduction  in  flutter  speed  associated  with  failure  of  the  outboard 
actuators  and  also  with  loss  of  stiffness  under  blow-back  and  stalled  actuator  conditions 
described  in  section  4.2  was  investigated  using  the  41  7,400-lb  configuration. 


4.4.3. 1 Failed  Outboard  Actuators 

Figure  18  shows  the  flutter  speeds  for  the  condition  where  the  outboard  actuators  are 
failed  and  the  inboard  actuators  are  rigid.  The  flutter  speed  boundary  and  modal  coupling 
for  the  failure  case  are  very  similar  to  the  nonfailure  case,  indicating  little  effect  from  the 
loss  of  the  outboard  links. 

When  the  actuators  and  associated  control  system  are  substituted  for  the  inboard  links 
(i.e.,  when  the  transfer  functions  of  figs.  8 and  ()  are  used),  the  flutter  boundaries  become  as 
shown  in  figure  19.  Two  significant  changes  in  the  flutter  characteristics  are  seen  by 
comparing  figure  19  with  the  rigid-actuator  configuration  of  figure  18.  First,  the  critical 
flutter  mode  is  now  the  first  body-bending  mode,  both  subsonically  and  supersonically. 
Second,  a considerable  decrease  in  flutter  speed  has  resulted,  again  most  significantly  in  the 
supersonic  regime.  There  remains,  however,  adequate  actuator  stiffness  to  prevent  stabilizer 
flutter  within  the  design  flight  placard,  which  is  for  the  failure  case. 

4.4.3. 2 Stalled  Actuators 

Figure  20  shows  the  flutter  boundary  when  the  stalled-actuator  transfer  function  (figs. 
10  and  II)  is  coupled  with  the  elastic  equations  of  motion.  The  stabilizer  flutter 
characteristics  are  similar  to  the  actuator-failure  condition  in  that  the  critical  flutter  mode 
originates  as  the  first  body-bending  mode  throughout  the  Mach  regime.  The  stalled-actuator 
transfer  function  exhibits  low  stiffness  in  the  flutter  frequency  region,  thus  resulting  in 
considerably  reduced  flutter  speeds.  Figure  20  shows  that  the  flutter  boundary  is,  however, 
above  the  design  flight  placard  for  all  Mach  numbers. 

4.4.3. 3 Blown-Back  Actuators 

The  blown-back  transfer  function  of  the  reference  3 tests  was  known  to  produce  an 
extremely  low  static  divergence  and  therefore  was  not  analyzed  in  detail  for  all  flight 
conditions.  The  proposed  blown-back  fix  was  an  attempt  to  remedy  the  divergence 
condition;  the  transfer  function  for  it  is  shown  in  figures  21  and  22.  As  can  be  seen,  the 
transfer  function  is  of  the  sixth  order  and  as  such  was  beyond  the  capabilities  of  the  flutter 
solution  computer  program  to  use  explicitly. 

An  approximation  to  the  proposed  blown-back  fix  transfer  function  was  developed 
with  a second-order  transfer  function  as  shown  in  figures  21  and  22.  Since  the  flutter 
frequencies  are  in  the  7.5-  to  9.5-rad/sec  range,  it  is  felt  that  the  approximate  second-order 
transfer  function  represents  the  proposed  sixth-order  transfer  function  satisfactorily  for  this 
blown-back  analysis. 

The  second-order  transfer  function  was  then  coupled  with  the  elastic  equations  of 
motion,  the  result  being  the  flutter  boundary  shown  in  figure  23.  Again,  the  critical  flutter 
mode  originated  as  the  first  body-bending  mode. 


The  proposed  blown-baek  fix  transfer  function  still  does  not  possess  adequate  stiffness 
to  prevent  flutter  in  the  supersonic  regime,  although  it  offers  a significant  improvement  over 
the  static  divergence  situation  produced  by  the  reference-3  transfer  function.  It  is  likely  that 


modification  of  the  system  to  decrease  the  aetuator  phase  angle  would  substantially  aid  in 
providing  flutter  clearance. 

4.5  NATURAL  MODES  AND  FREQUENCIES 

Table  2 summarizes  the  modal  frequencies  for  the  three  weight  conditions  analyzed. 
For  these  reference  mode  cases,  all  horizontal-stabilizer  actuators  are  represented  by  rigid 
links.  Also  noted  in  the  table  are  those  modes  retained  for  the  flutter  analyses,  i.c.,  those 
modes  that  contribute  significantly  to  stabilizer  flutter.  These  mode  shapes  are  plotted  and 
presented  in  reference  4. 

Table  3 presents  the  modal  frequencies  for  the  case  with  outboard  actuators  failed 
compared  with  the  case  with  all  four  actuators  connected.  Again,  for  these  reference  cases 
the  actuators  are  considered  to  be  rigid  links.  As  expected,  the  modal  frequencies  for  the 
two  configurations  are  nearly  identical  except  for  those  modes  involving  horizontal- 
stabilizer  motion.  Reference  4 contains  the  plotted  mode  shapes  for  the  two-actuator 
condition  on  the  lightweight  (417,400-lb)  airplane. 

4.6  FLUTTER  VELOCITIES  AND  FREQUENCIES 

Table  4 summarizes  the  flutter  speeds  and  frequencies  for  all  configurations  of  the 
actuator  that  were  analyzed.  Of  additional  interest  were  the  results  of  two  studies  made  at 
Mach  0.9  and  2.1  wherein  the  wing  aerodynamic  forces  were  deleted  from  the  flutter 
solution.  Generalized  inertia  and  stiffness  contributions  were  retained,  however.  Very  little 
change  in  the  horizontai-tail  pitch  flutter  mode  was  observed;  however,  the  body-bending 
flutter  mode  disappeared,  leading  to  the  conclusion  that  wing  aerodynamic  forces  are 
needed  for  a thorough  assessment  of  empennage  flutter  characteristics.  Subsequent  studies 
indicated  that  the  body-bending  flutter  mode  was  also  affected  significantly  by  control 
system  characteristics.  Nevertheless,  the  cantilevered  aft-body  analyses  can  provide 
significant  flutter  results  for  preliminary  design  assessment. 
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TABLE  3.-M0DAL  FREQUENCY  COMPARISON  FOR  SYMMETRIC  STABILIZER 

ACTUATOR  FAILURE 


Actuators 

Actuators 

Mode 

2 outbd 
2 inbd 

2 inbd 

Mode 

2 outbd 
2 inbd 

2 inbd 

1 

0 

0 

18 

9.786 

9 786 

2 

0 

0 

19 

10.67 

10.19 

3 

0 

0 

20 

11.38 

11.38 

4 

1.452 

1.452 

21 

11.52 

11.52 

5 

1.920 

1.920 

22 

12.89 

12,89 

6 

2.701 

2.701 

23 

14.04 

14.04 

7 

3.078 

3.078 

24 

14.33 

14.33 

8 

3.721 

3.721 

25 

14.66 

14.63 

9 

3.868 

3.868 

26 

14.74 

14.73 

10 

4.021 

4.021 

27 

14.95 

14.94 

11 

5.268 

5.267 

28 

15.26 

15.26 

12 

5.710 

5.707 

29 

15.54 

15.54 

13 

7.002 

6.989 

30 

15.81 

15.81 

14 

7.558 

7.536 

31 

16.63 

16.63 

15 

8.047 

7.975 

32 

1689 

16.89 

16 

9.203 

9.126 

33 

17.18 

17.17 

17 

9.478 

9.450 

34  1 

18.09 

18.09 

Free-free  airplane  modal  frequencies  (cps) 

Note;  rigid  actuators;  417,400-lb  gross  weight  airplane 
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ILIZER  SYMMETRIC  FLUTTER 


5.0  CANTILEVERED  AFT  BODY  ANALYSIS 


5.1  DESCRIPTION  OF  STRUCTURE 

A separate  study  was  performed  to  assess  the  effect  of  elevator-gearing  ratios  upon  the 
flutter  speeds  of  an  all-movable  stabilizer.  The  analytical  structural  model  of  the  fuselage 
used  for  these  analyses  consisted  of  an  idealization  of  the  SST  aft  body,  cantilevered  at  the 
wing  rear  spar.  The  structural  idealization  was  identical  to  the  corresponding  portion  of  the 
structure  used  for  the  free-free  airplane  analysis  described  in  section  4.0. 

Figure  24  shows  schematically  how  elevator  gearing  is  defined  and  accomplished.  The 
gearing  ratio  is  nearly  constant  for  small  stabilizer  rotation  angles  (i.e.,  when  6S«  0),  but 
deviates  from  this  value  near  the  design  limits  for  stabilizer  rot-.t'on  (i.e.,  when  6S«  ± 15°). 
Normal  modes  and  generalized  aerodynamic  forces  were  calculated,  however,  for  elevator 
gearing  at  fis  = 0. 

Three  gearing  ratios  were  examined  for  their  effect  on  stabilizer  flutter.  The  SST 
stabilizer  design  gearing  ratio  (tj)  of  1.77  was  analyzed,  as  well  as  V = 0 (no  gearing)  and 
tj=  2.5.  These  gearing  ratios  were  chosen  to  cover  the  range  of  values  considered  in  earlier 
designs  of  the  SST  stabilizer. 

Each  flutter  analysis  originated  with  the  generation  of  new  stiffness  matrices  for  the 
particular  gearing  ratio  being  examined.  For  the  non-geared  elevator,  this  preliminary  step 
consisted  of  removing  the  structural  members  representing  the  actuation  arm  and  link  and 
inserting  a localized  member  to  lock  the  elevator  to  the  main  stabilizer  structure.  Thus,  the 
stabilizer  and  elevator  rotated  as  a single  surface. 

A gearing  ratio  of  2.5  was  accomplished  by  modifying  the  geometry  of  the  actuation 
arm  and  link  structure.  Figure  25  shows  how  this  was  accomplished  and  shows  the  resulting 
change  in  elevator  rotation  in  comparison  with  the  configuration  with  the  gearing  ratio 
of  1.77. 

5.2  DESCRIPTION  OF  CONTROL  SYSTEM 

The  control  system  used  for  the  cantilevered  flutter  analyses  was  identical  to  that  used 
in  the  free-free  analyses  except  that  the  longitudinal  control  system  loops  (ECSS  and  USAS) 
were  not  included.  A block  diagram  of  the  control  system  is  shown  in  figure  26.  As  was 
done  in  the  free-free  analysis  of  section  4.0,  the  horizontal-stabilizer  actuators  were  initially 
assumed  to  be  rigid  links.  The  control  system  and  actuator  mode  equations  were  then 
merged  with  the  elastic  structure  equations  as  described  in  reference  4.  The  holding-actuator 
transfer  function  of  figures  7 and  8 was  used  when  the  control  system  effects  were  assessed. 
No  analyses  were  made  for  the  stalled  or  blown-back  actuator  conditions. 
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5.3  ANALYSIS  RESULTS 
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actuators  were  assumed  tor  the  stabilizer. 


When  (he neoto. con.ro, -system 

«*  ""  Varia,i0"  nU,'Cr 

speeds  being  on  the  order  ol  2.5  percent. 


i tivo  decreasing  the  rotational  constraining  stiffness 

|„  general,  it  was  observed  that  decree b . tt,„dej  IO  make  the  nutter 

'^rsn  w„h  .*»«.  ^ «*> ^h* 

increased  flutter  speed.. 


5.3.1  Natural  Modes  and  Frequencies 


Tabte  5 shows  the  modal . ST  ^ 

modes  were  used  in  all  nutter  soluti  . , er;iting  condition  (all  tour 

Ss^n;^rr^r"o„d,.,on  (outboard  actuators  faded,.  Referenced 

Shows  the  plotted  mode  shapes  for  each  analysis  configuratio  . 


5.3.2  Flutter  Speeds  and  Frequencies 


5. 3. 2.1  Rigid  Actuators 


C-  07  and  28  show  the  rigid  actuator  nutter  results  for  the  configurations  both 

with  S^t^n,,  re,  attach  = 0.9  and 

^riT^eiV  — * 


5. 3. 2.2  Control  System  Effects 


C,L  V/Uiiuwi  

including  the  con, to,  system  of 

condition  appears  to  amplify r”™w^  £ ring  r3tio  changes  is  still  small,  the  maximum 
speeds.  The  change  in  Hutte  p Inclusion  of  the  control  system  also  seems  to 

S2JS  abTeL°Ue  bend  in  the  effect  of  gearing  ra„o  on  Huttcr  speed;  f, gores  29  and  30 
show  that  increasing  the  gearing  ratio  tends  to  increase  the  nutter  speeds. 
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TABLE  5. -MODAL  FREQUENCIES  FOR  GEARED  ELEVATOR  FLUTTER  ANALYSES 


Modal  frequencies3  (cps) 


Gearing 

ratio 

529.9  kips 

690.3  kips 

Mode 

All 

act. 

oper. 

Outbd 

act. 

failed 

All 

act. 

oper. 

Outbd 

act. 

failed 

0 

1 

2.37 

2.37 

1.92 

1 92 

2 

10.21 

9.92 

7.85 

7.78 

3 

12.82 

12.48 

11.75 

11.29 

4 

14.76 

14  73 

14.16 

14.05 

5 

17.95 

17.85 

14  71 

14.67 

6 

20.43 

20.39 

19.00 

18.86 

7 

23.80 

23.69 

20.07 

20.07 

8 

24.58 

24.58 

21.72 

21.71 

1.77 

1 

2.369 

b 

1.92 

1.92 

2 

9.64 

1) 

7.78 

7.63 

3 

11.85 

b 

10.60 

10.05 

4 

14.82 

b 

13.59 

13.56 

5 

16.03 

b 

14.82 

14.79 

6 

20.05 

b 

16.67 

16.65 

7 

23.2 

b 

20.06 

20.06 

8 

24.57 

b 

22.10 

22.09 

2.5 

1 

2.37 

2.37 

1.92 

1.92 

2 

9.32 

8.70 

7.73 

7.53 

3 

11.60 

11.50 

10.16 

9.65 

4 

14.77 

14.74 

13.46 

13.44 

5 

15.73 

15.73 

14.82 

14.79 

6 

20.01 

19.99 

16.35 

16.35 

7 

23.15 

23.13 

20.06 

20.06 

8 

24.56 

24.56 

22.12 

1 22.12 

aAII  actuators  assumed  rigid,  and  fuselage  structure  cantilevered  at  wing  rear  spar 
^Modes  not  calculated 


6.0  CONCLUSIONS  AND  RECOMMENDATIONS 


The  following  summarize  the  conclusions  anil  recommendations  ol  this  study: 


1 ) The  analyses  performed  indicate  that  the  SSI  empennage  provides  adequate 
margins  ol  safety  for  flutter  for  all  flight  conditions  (including  hydraulic  actuator 
failure)  except  when  the  stabilizer  is  blown  back  due  to  high  aerodynamic 
loading. 

2)  The  analyses  further  indicate  that  the  horizontal-tail  flutter  speeds  are  relatively 
insensitive  to  elevatoi  gearing-ratio  changes. 


3j  Differences  in  flutter  speeds  obtained  from  the  analyses  conducted  under  this 
contract,  which  employed  updated  three-dimensional  unsteady  aerodynamics,  and 
from  corresponding  analyses  completed  previously  point  out  the  need  for 
experimental  confirmation  of  analytical  results  for  conligurations  of  this  sort.  It  is 
therefore  strongly  recommended  that  tests  of  the  I /7-scale  SST  transonic 
empennage  fluttei  model  be  conducted  to  provide  a basis  for  correlation  between 
analysis  and  test. 

4)  Although  the  linear  representation  of  the  hydraulic  system  stiffness  used  in  the 
analyses  is  considered  to  be  a satisfactory  approach  for  first-order  assessment  ol 
flutter  characteristics  of  systems  of  this  type,  further  work  to  develop  analytical 
methods  that  account  for  hydraulic  system  nonlinearities  is  recommended.  Sucu 
studies  should  employ  analog-digital  (hybrid)  computer  techniques  to  evaluate 
transient  actuator  behavior. 
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FIGURE  4.-HOR/ZONTAL  STABILIZER  STRUCTURAL  CONFIGURATION 


FIGURE  9.-H0RIZ0NTAL  STABILIZER  ACTUATOR  PHASE  ANGLE , 
HOLDING  CONDITION 


FIGURE  11. -HORIZONTAL  STABILIZER  ACTUATOR  PHASE  ANGLE, 
STALLED  CONDITION 
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FIGURE  12.—HORIZONTA L STABILIZER  ACTUATOR  STIFFNESS , 
BLOW  BACK  CONDITION 
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FIGURE  15 -SYMMETRIC  HORIZONTAL  STABILIZER  FLUTTER  BOUNDARIES, 
ACTUA  TORS  HOLDING 
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FIGURE  16.— AIRPLANE  GROSS  WEIGHT  EFFECTS  ON  SYMMETRIC  HORIZONTAL 
STABILIZER  FLUTTER 
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FIGURE  17-STABILIZER  FLUTTER  SPEED  VARIATION  WITH  TANK  8A  FUEL 
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FIGURE  20.-SYMMETRIC  HORIZONTAL  STABILIZER  FLUTTER  BOUNDARY 
ACTUATORS  STALLED 


FIGURE  22.-TRANSFER  FUNCTION  APPROXIMA  TION,  ACTUA  TORS 
BLOWN  BACK 


Flutter  speed  (kt  EAS) 


FIGURE 23.-SYMMETRIC  HORIZONTAL  STABILIZER  FLUTTER  BOUNDARY, 
ACTUA  TORS  BLOWN  BACK 
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FIGURE  25-GEOMETRY  CHANGES  TO  PRODUCE  GEARING  RATIO  CHANGE 


Cantilevered 
structure 
equations  of 
motion 
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Transfer  function  varies  with 
actuator  condition.  Single,  holdina 
actuator  shown  acting  at  a distance 
of  33.1  in.  about  hinge 
< <5  s = 0) 


j = actuator  linear  displacement 
A divided  by  torque  arm  <33.1  in.) 

q = actuator  force  multiplied  by 
A torque  arm  (33.1  in.) 


FIGURE  26.-CONTROL  SYSTEM  USED  FOR  GEARED  ELEVATOR 
FLUTTER  ANALYSES 
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FIGURE  2B.-EFFECT  OF  ELEVA  TOR  GEARING  ON  STABILIZER  FLUTTER, 
OUTBOARD  ACTUATORS  FAILED,  INBOARD  RIGID 
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APPENDIX  A 


EQUATIONS  OF  MOTION  DEVELOPMENT 


The  theoretical  development  of  the  equations  of  motion  for  an  elastic  airplane  struc- 
ture including  a stiffness  coupled  actuator  mode  is  given  as  follows: 

Assume  a symmetric  analysis,  for  which  the  Lagrange  equations  can  be  written  as 

= Z Zjjp  <A1) 

t\3q  / j 


A I 

dt' 


where: 

T 

Z; 


= kinetic  energy 

= sum  of  structural  stiffness  force  plus  external  aerodynamic  forces 
= total  displacement  of  panel  i 
= generalized  coordinate 


Writing  out  the  above  expressions. 


* ! i-n 


(A2) 


Zi  " " 2 kij  zjE  2 p ? aij  aj 
J J 


where 


a,- 


- fi  ^ 

" U"  \9x /j 


ay  = aerodynamic  influence  coefficients 

ky  = stiffness  matrix  for  structure 

Zj^  = elastic  component  of  total  displacement 

In  terms  of  generalized  coordinates,  we  can  write 

z;  = 0jq 


(A3) 


(A4) 


(A5) 


Assume,  initially,  that  we  are  considering  only  the  elastic  modes  of  the  structure,  i.e\,  no 
actuator  mode  is  yet  included.  Therefore, 


Zj  = Zj 


(A6) 
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Introducing  the  superscript  F in  ail  terms  to  denote  the  elastic  components  only,  we  ct 
write  the  Lagrange  equations  as: 


can 


d_ 

dt 


i E 2) 


3q 


TT 


9z 


= -?  ? kijZF— i ipU2  22a.. 

' J 1J  J 3qt  2 " j j ij 

Eliminating  the  subscripts  and  writing  in  matrix  form, 

jzE}=  [0E)  jqEj 


U \ 9 X / : 


1 3z  E 


9q 


( A7 ) 


and 


the  above  equation  becomes. 


T 


10E1  [ni]  [0E]  jqEJ  = ~[0E]  [k]  [0E]{qE)--ipU[0E]T[a)  |0E]{qE} 


+ypU2  (0E)  (a)  (0E)  |qEJ 

Equation  (A8)  is  recognized  as  the  usual  equation  of  motion  for  an  aeroelastic  structure. 

Let  us  include  now  an  additional  actuator  mode.  The  effect  of  this  mode  is  to  create 
additional  deflections,  z A,  so  that 


(A8) 


7 -4  Z E 
i L j 


but, 


Zj  = z f + z A 

In  terms  of  an  actuator  generalized  coordinate,  5A,  we  can  write 


zt  = 5A 


(A9) 


(A  10) 


(All) 


As  stated  in  (A9)  above, 


jzj^  I0E)  |qEJ 


(A  1 2 ) 


m 


but  we  can  say  that 


\'\=  I0hl {q} 


< A 1 3 ) 


where  {q} is  now  a new  set  of  generalized  coordinates  that  relate  the  elastic  structure  mode 
shapes  to  the  total  displacements  of  the  structure  nodes. 


Extend  the  bar  notation  to  all  terms  to  denote  the  inclusion  of  the  actuator  mode,  so 
that  (A  13)  becomes 


fz|=  ly>L  l{<)} 


(A  1 4) 


Lagrange's  equations  now  appear  as, 


_d_ 

dt 


3(f  2mi?i2) 


3(t 


= _ 2 
i 


, d(7i  ~ ziA) 


- 2 2 -Ipu-a: 


-(I). 


3Z: 


bi  <A|5» 


fkij('rzjA)  a,r  , j 

Since  §A  is  now  another  generalized  coordinate,  another  equation  o(  motion  is  produced: 

a(i 


jL 

dt 


"aT 


= -2  2 k (z--z  A\^i~ZiA) 
i j K'J\J  7)  ) 35a 


- £ ^pU2a- 

i j *• 


U 


U Vax/: 


a?.: 


(A  1 6) 


'jj  aSA 

111  matrix  form,  equations  (A  15)  and  (A16)  may  be  combined  as 


<?A 


0* 


0 


[ml 


&E|°]fe-j+^u  -6ET'31  PH  j 


L(l_*  + 


jL 


-V 


ki  pi-4a1 


[k] 


1 , p 

-~P  u- 


(A  I 7) 


Nott  that  the  coupling  between  the  elastic  modes  and  the  actuator  mode  occurs  in  the  generalized 
stiffness  matrix,  thus  giving  rise  to  the  term  “stiffness  coupled  formulation.” 


Equation  (A1 7)  can  be  written  in  terms  of  unsteady  aerodynamics  for  use  in  flutter  studies. 


I Mi° 
I 0|0 

i 


°A 


JL jr.l^J T Ik]  I0AJ 

^ ! T 


L-  [0A1  Ik]  [0E]|  [0A1  [k]  10A) 


±. 


Plain 


A|0 


Oj  0 

I 


L 9lJ 


0 


■ 'Oa(  IAI8) 
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where: 


M = 
K = 
A = 


[fi] 

T [ml  [ 

0E' 

T[k]  [ 

r f1 

T f 

KJ 

[a]  |_ 

(A18) 


U 


0E-0'E 


We  can  make  Qa  a generalized  coordinate  and  introduce  structural  damping  to  the 
equations  of  motion  in  the  usual  manner.  Equation  (A18)  then  becomes 


Oi  0:0 


K(_l  + jg)|_K  1 2j  _o_ 

Kt i ! Kii 


*|  / ’j 

( \ 

ilI 

2 

9 ( 

7 l 

L‘a  I 

a;oio 

oloTo 

,:d 

JIoAl 

L j 

0A) 

= {o} 


(A19) 


where: 

*12 

K' 


“21 


K 


22 


- [^E] T i k i M 

-[0A]T[k]  [d>E] 

[0A]T  [kl  [</>A] 


The  control  system  feedback  transfer  function  relating  the  actuator  torque  and 
displacement  5^  is  generally  of  the  form: 


where: 


T(s)  = 


_ Qa  ^2  (s  + ) 


S + cj' 


C2,  <0  j,  a>2  = parameters  of  the  system  design 

S = Laplace  transform  operator 

In  order  to  merge  the  control  system  equation  with  the  equations  of  motion,  it  is 
necessary  to  transform  the  equations  of  motion,  (A19),  giving 


, 1 1 

M^O^O 

“0;  o;o~ 

s2  + 

K(1  +ig)j  k12:  0 
^21  1^22  [-••0 


0 ' 

aLoio 

To" 

+ r 

ojojo 

= [0] 


(A20) 
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The  control  system  transfer  function  can  now  be  merged  with  the  equations  ot  motion 
to  form. 


MjO'O 
0 |0|  0 
o”io!o 


o;  o! 

0 

K(l  +ig)|  Kpio 

9 

a!  o ; o 

i | 

s2  + 

*~i — {■ 
O'  0. 

' -4 

0 

s + 

^21  ! ^22  ! 

+ T 

o ! o j o 

_°!‘c2! 

-1.0 

0 j i j -co-) 

_0|0|0_ 

= 0 


(A21 ) 


Simplifying  the  notation, 

2 

[ B]  S2  + fCl  S+  [D]  + ^ |E]  =0  (A22) 


This  is  the  equation  to  be  solved  for  flutter  instabilities.  Note,  however,  that  both  S and 
o>  are  unknown  complex  variables.  If  we  restrict  our  examination  of  the  solution  to  the 
vicinity  of  the  imaginary  axis  ( a = 0 where  S = a + ico)  we  can  write  equations  (A22)  as 


[B]  --pin] 


S2  + |C]  S + [D|  =0 


(A23) 


which  is  then  rooted  for  S.  From  these  roots,  flutter  velocity,  frequency,  and  damping  values 
can  be  obtained. 


f<-  = flutter  frequency  =rr-  ABS  (ROOT) 

f = damping  = RKAL  (ROOT)/27rf|- 

g = -2.0  (f)  (for  small  damping) 

27Tbrff 

Vf  = flutter  velocity  = — £ — 


where  k is  the  reduced  frequency  .ssociated  with  the  airforce  matrix  [H] . 
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